基于高斯烟羽模型的船舶尾气扩散研究

为有效对港区大气污染进行治理、分析船舶尾气,本文详细介绍了一种基于高斯烟羽模型,通过MATLAB模拟仿真模型,其包括实验仿真过程、技术原理及理论模型对船舶尾气扩散进行的研究。该模型是在传统的高斯烟羽模型的基础上,通过对实源像源进行加权选择输入参数;通过矢量合成确定了气体扩散的方向,利用合成后的"风速"进行计算仿真,有效模拟了船舶尾气在港区或者海洋环境中的气体扩散模型。其模型简单且可以有效模拟船舶尾气扩散。并且进一步对后续模型的精确优化进行分析。     

基于分布式控制力矩陀螺的水下航行器轨迹跟踪控制

基于控制力矩陀螺群(CMGs)的水下航行器具有低速或零速机动的能力。采用基于分布式CMGs的水下航行器方案,并研究其水平面的轨迹跟踪控制问题。通过全局微分同胚变换将非完全对称的动力学模型解耦成标准欠驱动控制模型,并根据简化的模型构建其轨迹跟踪的误差动力学模型,将轨迹跟踪控制问题转化为误差模型镇定问题。基于一种分流神经元模型和反步法设计了系统的轨迹跟踪控制律,该控制器不需要对任何虚拟控制输入进行求导计算,且能确保跟踪误差的最终一致有界性。仿真结果表明该控制器能够实现在不依赖动力学参数先验知识的情况下对光滑轨迹的有效跟踪。     

Ambipolar diffusion in the solar atmosphere

The process of non-linear ambipolar diffusion in the region overlying the solar surface can be an effective mechanism for producing sharp magnetic structures and current sheets. These may be the sites responsible for the occurrence of connectivity of magnetic field lines, and the subsequent explosive input of energy for heating of some of the features in the atmosphere of the Sun..     

MONITORING THE DIFFUSION OF SUSPENDED SEDIMENTS AND STABILITY OF TIDAL RADIAL SAND RIDGES AREA USING MULTI-SATELLITE REMOTE SENSING DATA

Nine Landsat TM tapes and images and MSS images, 10 NOAA tapes and images and 1 SAR image from 1973 to 1997 were used to analyse the diffusion of suspended sediments and the change of tidal radial sand ridges in the northern part of the Changjiang River delta, the South Yellow Sea. The results showed that the diffusion of suspended sediments was controlled by the tide, net current, and submarine topography in this area. The distribution of suspended sediments had close relationship with thesubmarine topography. The old Huanghe River delta and the Changjiang River comprise the main sediment supply for the formation of radial sand ridges, whose evolution can be divided into three stages since the Huangbe River changed its course and flowed northward into the Bohai Sea.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Asymmetrically zoned reaction rims: assessment of grain boundary diffusivities and growth rates related to natural diffusion-controlled mineral reactions

This study explores garnet coronas around hedenbergite, which were formed by the reaction plagioclase + hedenbergite→garnet + quartz, to derive information about diffusion paths that allowed for material redistribution during reaction progress. Whereas quartz forms disconnected single grains along the garnet/hedenbergite boundaries, garnet forms ～20‐μm‐wide continuous polycrystalline rims along former plagioclase/hedenbergite phase boundaries. Individual garnet crystals are separated by low‐angle grain boundaries, which commonly form a direct link between the reaction interfaces of the plagioclase|garnet|hedenbergite succession. Compositional variations in garnet involve: (i) an overall asymmetric compositional zoning in Ca, Fe²⁺, Fe³⁺ and Al across the garnet layer; and (ii) micron‐scale compositional variations in the near‐grain boundary regions and along plagioclase/garnet phase boundaries. These compositional variations formed during garnet rim growth. Thereby, transfer of the chemical components occurred by a combination of fast‐path diffusion along grain boundaries within the garnet rim, slow diffusion through the interior of the garnet grains, and by fast diffusion along the garnet/plagioclase and the garnet/hedenbergite phase boundaries. Numerical simulation indicates that diffusion of Ca, Al and Fe²⁺ occurred about three to four, four and six to seven orders of magnitude faster along the grain boundaries than through the interior of the garnet grains. Fast‐path diffusion along grain boundaries contributed substantially to the bulk material transfer across the growing garnet rim. Despite the contribution of fast‐path diffusion, bulk diffusion through the garnet rim was too slow to allow for chemical equilibration of the phases involved in garnet rim formation even on a micrometre scale. Based on published garnet volume diffusion data the growth interval of a 20‐μm‐wide garnet rim is estimated at ～10³–10⁴ years at the inferred reaction conditions of 760 ± 50 °C at 7.6 kbar. Using the same parameterization of the growth law, 100‐μm‐ and 1‐mm‐thick garnet rims would grow within 10⁵–10⁶ and 10⁶–10⁷ years respectively.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Possible Density Change of the Origin of North Pacific Intermediate Water Due to Double Diffusion in the Mixed Water Region

In this study we test Talley's hypothesis that Oyashio winter mixed-layer water (26.5–26.6σ _θ) increases its density to produce the North Pacific Intermediate Water (NPIW) salinity minimum (26.7– 26.8σ_θ) in the Mixed Water Region, assuming a combination of cabbeling and double diffusion. The possible density change of Oyashio winter mixed-layer water is discussed using an instantaneous ratio of the change of temperature and salinity along any particular intrusion (R _l ). We estimate the range of R _l ^DD required to convert Oyashio winter mixed-layer water to the NPIW salinity minimum due to double diffusion, and then assume double-diffusive intrusions as this conversion mechanism. A double-diffusive intrusion model is used to estimate R _l ^DD in a situation where salt fingering dominates vertical mixing, as well as to determine whether Oyashio winter mixed-layer water can become the NPIW salinity minimum. Possible density changes are estimated from the model R _l ^DD by assuming the amount of density change due to cabbeling. From these results, we conclude that Oyashio winter mixed-layer water contributes to a freshening of the lighter layer of the NPIW salinity minimum (around 26.70σ_θ) in the MWR.     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.